Nature of the Quantum Phase Transition in Clean, Itinerant Heisenberg Ferromagnets 



T.R.Kirkpatrick 

Institute for Physical Science and Technology, and Department of Physics 
University of Maryland, College Park, MD 20742 



m 
o 
o 

(N 



a 

o 



> 
o 

ON 

o 

(N 

o 



X3 
O 

o 



X 



D.Belitz 

Department of Physics and Materials Science Institute, University of Oregon, Eugene, OR 97403 

(Dated: February 1, 2008) 

A comprehensive theory of the quantum phase transition in clean, itinerant Heisenberg ferromag- 
nets is presented. It is shown that the standard mean-field description of the transition is invalid 
in spatial dimensions d < 3 due to the existence of soft particle-hole excitations that couple to the 
order parameter fluctuations and lead to an upper critical dimension d^ = 3. A generalized mean- 
field theory that takes these additional modes into account predicts a fluctuation-induced first-order 
transition. In a certain parameter regime, this first-order transition in turn is unstable with re- 
spect to a fluctuation-induced second-order transition. The quantum ferromagnetic transition may 
thus be either of first or of second-order, in agreement with experimental observations. A detailed 
discussion is given of the stability of the first-order transition, and of the critical behavior at the 
fluctuation-induced second-order transition. In d = 3, the latter is mean field-like with logarithmic 
corrections to scaling, and in d < 3 it can be controlled by means of a 3 — e expansion. 

PACS numbers: 75.40.Cx: 75.40.Gb; 64.60Ak 



I. INTRODUCTION 

The study of quantum phase transitions (QPTs) is cur- 
rently an important and very active field of research in 
condensed matter physics, see, e.g., Ref. ? . Although, 
strictly speaking, these transitions occur only at zero 
temperature [T — 0), they are important for understand- 
ing the behavior of many systems at low, but routinely 
accessible, temperatures. Understanding QPTs is also 
important for gaining insight into the possible phases of 
systems at zero temperature. Indeed, QPTs are thought 
to be relevant for understanding phenomena as diverse as 
high-Tc superconductivity, the quantum Hall effects, var- 
ious magnetic phenomena in both metallic and insulat- 
ing systems, the transport properties of doped semicon- 
ductors, and superconductor-metal and superconductor- 
insulator transitions, see, e.g., Refs. 1111. 

Perhaps the most obvious, and one might naively think 
almost trivial, QPT is the ferromagnetic transition that 
takes place in a clean^ itinerant electron system as the 
exchange coupling is varied at zero temperature. Indeed, 
this was one of the first QPTs to be considered, see Ref. 
? and references therein. The traditional arguments 
and results for this QPT can be paraphrased as follows. 
Let M be the order parameter vector, i.e., the mag- 
netization, with components (i = x,y,z). Landau 
theory,^ which as a mean-field description is suitable for 
both thermal and quantum phase transitions, says that 
the free energy (at T > 0), or energy (at T = 0), as a 
function of M , for small magnetization and small mag- 
netic field h is of the form 



h-M 



(1.1) 



Here Fq is a background contribution that reflects the 
degrees of freedom other than the order parameter, t 



turns out to be the distance from the mean-field transi- 
tion, i.e., the transition takes place at i = 0, and u is 
a constant that is assumed to be positive. Minimizing 
Eq. Hl.lfl with respect to M leads to the conclusion that 
at zero external magnetic field the magnetic transition is 
continuous with mean-field or Landau critical exponents. 
Fluctuations invalidate Landau theory in sufficiently low 
dimensions, while the mean- field critical behavior is exact 
in dimensions d larger then an upper critical dimension 
+ . For the thermal phase transition in a Heisenberg 
ferromagnet, it is well established that (i+ = 4. For the 
corresponding QPT, it was argued that — A — z, with 
z the dynamical scaling exponent.'^ This reduction of 
the upper critical dimension is a result of the coupling 
between statics and dynamics in quantum statistical me- 
chanics, which leads to an effective dimension for fluctu- 
ations given by deff = d + z. Mean- field theory suggests 
2 = 3 for the quantum Heisenberg transition of clean 
itinerant electrons, so the conclusion was that this QPT 
should have a d^ — 1, resulting in mean field- like critical 
behavior for both thin films and bulk systems. From a 
theoretical statistical mechanics point of view, the itiner- 
ant quantum ferromagnetic transition therefore did not 
appear to be very interesting. 

This conclusion was recently challenged by what 
amounts to a generalized mean-field description of the 
transition.'" The basic physical argument, which is gen- 
eral and applies to other phase transitions as well, is as 
follows. In the disordered phase, Fq contains contribu- 
tions from fermionic soft modes, viz., particle- hole ex- 
citations. Some of these acquire a mass in the ordered 
phase, which decreases the contribution of these fluctu- 
ations to the free energy, and thus leads to a negative 
term in the free energy function that has a nonanalytic 
dependence on the order parameter If this mode-mode 
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coupling effect, which is neglected in the usual Landau 
or mean-field theory, is strong enough, it clearly can lead 
to a modification of phase transition predicted by Eq. 

Ref. ? showed that in the case of an itinerant 
fcrromagnet, the soft modes that couple most strongly to 
the order parameter, viz., spin-triplet particle-hole exci- 
tations, do indeed develop a mass in the ordered phase 
and lead to a Landau energy function that has the form, 
in d = 3, 

F = Fo+tM^ +vM'^\iiM^ +uM^ -h- M . (1.2) 

The presence of the M'^ In Af ^ term, compared to Eq. 

changes the nature of the transition from a con- 
tinuous one to a discontinuous one. The same is true 
for dimensions 1 < d < 3, due to a similar nonanalytic 
term in the Landau function. The fact that the nature of 
the phase transition in d = 3 changes qualitatively upon 
improving on Landau theory is not consistent with the 
traditional notion of a 3-d system being above its upper 
critical dimension = 1. In contrast to the traditional 
prediction of a continuous transition with mean-field ex- 
ponents, the generalized mean-field theory predicts the 
transition to always be of first order provided that d < 3. 

Experimentally, the situation is seemingly inconclu- 
sive. In some ferromagnets with low Curie temperatures, 
where the quantum phase transition can be triggered by 
hydrostatic pressure, or composition, the observed tran- 
sition is of first-order, in agreement with the general- 
ized mean-field theory. This is the case, for instance, in 
MnSi,' and in UGe2.'^ In others, for instance, ZrZn2,'^ 
and NixPdi_x,'^ however, the transition is observed to 
be continuous. Moreover, the critical behavior observed 
in NixPdi_x is in good agreement with mean-field expo- 
nents. This is surprising, given the above conclusion that 
mean- field theory cannot be correct in d = 3. 

In this paper we provide new insights into this QPT, 
and additional understanding of the discontinuous tran- 
sition that results from the generalized mean-field theory, 
as well as of the stability of the latter. In the general the- 
ory of phase transitions, transitions that are predicted to 
be continuous by Landau theory but are in fact discon- 
tinuous are called fluctuation-induced first-order phase 
transitions.'^ We will show that the first-order transition 
in itinerant quantum ferromagnets can indeed be under- 
stood as being fluctuation-induced. The novel feature is 
that although the order parameter fluctuations are above 
their upper critical dimension, in a well-deflned sense the 
soft fermion fluctuations are not, and it is the effect of 
these fluctuations that drives the transition first order for 
d < 3. 

Their are many similarities between the fiuctuation- 
induced first-order phase transition discussed here, and 
the thermal first-order transition that occurs in con- 
ventional superconductors, or the nematic-to-smectic-A 
transition in liquid crystal systems.'^ ' In all of these sys- 
tems there are soft or massless excitations (in supercon- 
ductors, these are the gauge or vector potential fluctua- 
tions; in liquid crystals, the director fluctuations; while in 



the electron system considered here, fermionic particle- 
hole fluctuations) that couple to the order-parameter 
fluctuations and become massive in the broken-symmetry 
phase. Because of the latter property, the fluctuation 
contribution to the free energy decreases, which ulti- 
mately leads to a fluctuation-induced flrst-order transi- 
tion. If these fluctuations are integrated out in some ap- 
proximation, then a nonanalytic Landau-like theory can 
be derived which predicts a discontinuous phase transi- 
tion. The modified mean-field theory for the magnetic 
transition mentioned above results from such a proce- 
dure. For the superconductor and liquid crystal transi- 
tions, a similar technique was used initially.'" Later, a 
renormalization-group description of fiuctuation-induced 
thermal phase transitions was developed.'^ Part of our 
goal here is to do the same for the quantum ferromag- 
netic transition. 

In addition, we perform a renormalization-group anal- 
ysis of the stability of the first-order transition predicted 
by the generalized mean-field theory. It turns out that 
the first-order transition is stable if it occurs at a suffi- 
ciently large value of the paramter t in Eq. H1.2() . How- 
ever, if it occurs at small values of t, then the first-order 
transition can in turn become unstable with respect to 
fluctuations. The final result in that case is a second- 
order transition that is induced by fiuctuations, by a 
mechanism that is similar to the one discussed in the con- 
text of classical Potts models by Fucito and Parisi.'^ This 
second-order transition is distinct from Hertz's mean- 
field transition and belongs to a different universality 
class. Depending on microscopic parameter values, the 
ferromagnetic QPT in itinerant electron systems can thus 
be either of first order, or of second order, in agreement 
with the experimental observations mentioned at the be- 
ginning of this section. Moreover, the critical behavior in 
the continuous case in d = 3 is mean field-like with loga- 
rithmic corrections to scaling. Within the current experi- 
mental accuracy, this is indistinguishable from mean-field 
exponents, again in agreement with the experimental ob- 
servations. 

The outline of this paper is as follows. In Section ^ 
we give a basic field theory describing coupled magne- 
tization fluctuations and soft fermionic degrees of free- 
dom. In Section lTTTl we derive and discuss the generalized 
mean-field theory that results from integrating out the 
fermionic modes and that predicts a first-order transition. 
In Sec. IIVI we perform a one-loop renormalization-group 
analysis of the field theory. We show that the renormal- 
ized quartic coupling constant can become negative at 
large scales, leading to the first-order phase transition 
described by the generalized mean-field theory. How- 
ever, in a certain parameter regime fluctuations keep the 
quartic coefficient positive, which results in a continuous 
phase transition. In Section we further discuss our re- 
sults, compare them with previous work, and comment 
on the experimental situation. Parts of the results pre- 
sented here have been previously announced in two short 
publications.'^ ^ 
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II. THE COUPLED FIELD THEORY 

Recently we have derived and discussed a local field 
theory describing the quantum ferromagnetic transition 
in disordered itinerant electron systems.^ ^ This field 
theory describes the coupling between the soft or mass- 
less fermionic degrees of freedom (which in a disordered 
electron system are diffusive, i.e., the frequency is a 
quadratic function of the wavenumber) and the magne- 
tization fluctuations. Here we give the analogous theory 
for clean electronic system. In this case the fermionic 
soft modes have a linear dispersion, i.e., the frequency is 
a linear function of the wavenumber. There are numerous 
ways to construct field theories that describe these soft 
modes; here we choose the method developed in Refs. ? 



A. Soft Modes 

A generalized Landau-Ginzburg- Wilson (LGW) theory 
that keeps all of the soft modes in the problem will be 
described in terms of an action A that depends on a field 
M describing the fluctuating magnetization, and on a 
field q describing the soft fermionic two-particle modes. ^ 
All other modes we imagine have been integrated out 
in order to arrive at the effective theory. The partition 
function can then be written in terms of a functional 
integral with respect to M and q, 



I D[M,q] 



(2.1a) 



The action will consist of a part that depends only on the 
magnetization, a part that depends only on the fermionic 
degrees of freedom, and a coupling between the two. 



A[M, q] ^AM + Aq + A 



M,q 



(2.1b) 



The various pieces of the action in Eq. (|2.1|) can be de- 
rived starting from a microscopic fermionic action, or 
more generally written down on the basis of symmetry 
arguments^ Here we choose the latter approach, with oc- 
casional references to the miscroscopic theory as a check. 
For a complete derivation from a microscopic action, the 
techniques of Refs. ? and ? can be used. 

Am is just a static, local, LGW functional for the mag- 
netization fluctuations. It is local because no massless 
modes that couple to the magnetization have been inte- 
grated out, and it can be chosen static because the rel- 
evant (in the long-wavelength, low-frequency limit) dy- 
namical part will be shown to be provided by the coupling 
to the q fluctuations, see Eq. H2.9c|l below. Neglecting 
terms that are irrelevant for our purposes. Am is given 
by 



Am = - j dx M{x) [t - aV^] M{x) -u j dx M^{x) 



Here x = {x,t) comprises the real space position x 
and the imaginary time r. J dx = J dx dr with 
P — l/ksT, where T is the temperature, t is the di- 
mensionless distance from the bare LGW critical point, 
and a and u are positive constants. The physical magne- 
tization m is proportional to the expectation value of the 
field M. For later reference we also define a temporal 
Fourier transform of the field M by 



Mn{x) = Vt dr e'""^ M{x) , 
Jo 



(2.3) 



with Qn — 27rTn a bosonic Matsubara frequency. 

The soft fermion field q originates from the composite 
fermion variables' 



/ 



-V'lt'02T 

-_V'uV'2T 

-01iV'2T 

-V-'iT'V^ZT 



-1pl-\1p2l 

-V'uV'2i 

--01|-02i 



i/'itV'2T 

V'liV'2T 

-_V'uV'2T 

'01tV'2T 

(2.4a) 

Here the ^ and ip ^re the fermionic, i.e., Grassmann- 
valued, fields that provide the basic description of the 
electrons,'^ and all fields are understood to be taken at 
position x. The indices 1,2, etc. denote the dependence 
of the Grassmann fields on fermionic Matsubara frequen- 
cies ujni — 27rT(ni -I- 1/2), etc., and the arrows denote 
the spin projection. It is convenient to expand the 4x4 
matrix in Eq. (|2.4a(l in a spin-quaternion basis,' 



Qi2{x) 



r,i—0 



{Tr «> Si)lQj^2i^) 



(2.4b) 



with To = So = I2 the 2x2 unit matrix, and Tj = —Sj = 
—icTj [j = 1,2,3), with cri,2,3 the Pauli matrices. In this 
basis, z = and i = 1,2,3 describe the spin-singlet and 
spin-triplet degrees of freedom, respectively. The r = 0, 3 
components correspond to the particle- hole channel (i.e., 
products V'V' or V'"'/')? while r = 1, 2 describe the particle- 
particle channel (i.e., products tp^ or V'V')- For our pur- 
poses the latter are not important, and we therefore drop 
the r = 1,2 from the spin-quaternion basis. In terms of 
the remaining fields, the spin density can be expressed as 



{x) CTab V'm+n.b(a;) 

m ab 

= ^ (\/^)''tr [{Tr®Si)Q„,,,n+n{x)] , 

(i = l,2,3) (2.4c) 



m r=0,3 



(2.2) 



The matrix elements of Q are bilinear in the fermion 
fields, so Q-Q correlation functions describe two-fermion 
excitations. In a Fermi liquid, the Q-fiuctuations are 
massive and soft, respectively, depending on whether the 
two frequencies carried by the Q field have the same sign, 
or opposite signs, respectively. We therefore separate 
the Q fluctuations into massless modes, (712, and mas- 
sive modes, P12, by splitting the matrix Q into blocks in 
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frequency space, 

Qnmix) = Q{nm) Pnrnix) + 6 (n) 9 (-m) (7„,„ (a;) 

+ei-n)&im)ql^ix) . (2.4d) 

In what follows, we will incorporate the frequency con- 
straints expressed by the step functions into the fields P 
and q, respectively. That is, the frequency indices of P 
must always have the same sign, and those of q and q^ 
must always have opposite signs. 

Finally, we define spatial Fourier transforms by 

M„(fc) = ^ J dx Mn{x) , (2.5) 
and analogously for the fields q and q^ . 

B. Soft-mode field theory 

The massive modes P can be formally integrated out 
to obtain an effective action for the soft modes, 
This can be done perturbativcly, as the P-dependent part 
of the action takes the form of a stable Gaussian (i.e., 
quadratic in P) piece, and terms of higher order in P as 
well as terms coupling P and q, and P and M, respec- 
tively, starting with bilinear coupling terms. As can be 
seen from Eq. H2.4a|l . the q are particle- hole excitations, 
which in a clean electron system have a linear disper- 
sion relation, i.e., the frequency scales linearly with the 
wavenumber. The Gaussian part of the fcrmionic action 
will therefore have the form 

1,2,3,4 r,i 

(2.6a) 

The vertex function F*^^^ is most easily written in mo- 
mentum space, 

Tf2^34(fe) = -513 62i rf2°'(fc) + <5i_3,2-4 ^^0 27:TGK, , 

(2.6b) 

with 

T^i2°\k) = \k\ + GHni^2 ■ (2.6c) 

Here G and H are model-dependent coefficientsi^ If one 
derives Eq. Ij2.6b|) from the microscopic model of Ref. ? 
, one finds G = tt'^Np/2vf and H = l/irNp, with vp the 
Fermi velocity, and the density of states per spin at 
the Fermi surface. More generally, however, G and H 
will be arbitrary coefficients with the appropriate dimen- 
sions. Ks is a spin-singlet interaction amplitude that we 
include in our Gaussian theory in a RPA-type fashion. 
Inverting F*^^) shows that its presence does not change 
the frequency-momentum structure of the g-propagator, 
see Eqs. (|2.10|l below. There is no spin-triplet interac- 
tion in the bare action since its effects are included in 
Am- In a formal derivation from a microscopic action, 
this can be achieved by means of a Hubbard-Stratonovich 



decoupling of the spin-triplet interaction, with M the 
Hubbard-Stratonovich field. ^ ^ However, as long as Ks 
is nonzero, it generates a spin-triplet interaction in per- 
turbation theory. This has important consequences, see 
Sec, nil A1 below, and it is the reason we include K^. 

The part of the action coupling M and q originates 
from a term Am-q that couples M and Q. Such a term 
must be present since in the presence of a magnetization 
the fermionic spin density will couple linearly to it. Using 
Eq. I|2.4c|l . we thus obtain 

3 

Am-q = 2ciVt da;^^M;(a;) 

n 1=1 

X (-l)''/^^tr[(T^«)Sj)Q™,„+„(a;)] , 

r— 0,3 m 

(2.7a) 

with a model-dependent coefficient ci. In a technical 
derivation from a model with a point-like spin-triplet in- 
teraction amplitude Kt, this term also is produced by a 
Hubbard-Stratonovich decoupling of the spin-triplet in- 
teraction term,'^ and ci = 7TKt/2. Defining a sym- 
metrized magnetization field by 

bi2{x)=Y,irr®s,)lb^^ix) , (2.7b) 

with components 

lb,,ix) = (-)^/2^J„,„,_„jM;(a;) 

n 

+ i-Y+HPL^{x)] , (2.7c) 
allows to rewrite Eq. H2.7a|l in a more compact form, 

Am-q = ciVr J dx tr {b{x) Q{x)) . (2.7d) 

Using Eq. (|2.4d|l in Eq. H2.7a|l or (|2.7d|l . and integrating 
out the massive P-fluctuations, obviously leads to a series 
of terms coupling M and q, M and g^, etc. We thus 
obtain Am.q in form of a series 

-4-M,g = AM-q + AM-q2 + ■■ ■ (2.8a) 

The first term in this series is obtained by just replacing 
Q by g in Eq. (|T7ajl . 

AM-q = ciPi/^y" dx tr {b{x) qix)) 

3 

= 8ciTl/2^ / dxY,Y.rbl2i^)r<ll2i^) • 

12 r i=l 

(2.8b) 

The next term in this expansion must have the overall 
structure 

AM-q^ / dx tr (b{x) q{x) q\x)) 
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The details require information about the structure of the 
massive modes that were integrated out in going from Q 
to q. From the derivation of the nonhnear sigma model 
that results in the disordered case if one integrates out 
the massive P-fluctuations in tree approximation'^ it 
is known that the resulting effective fermion matrix field 
is traceless, i.e., {q'^)nm in the above expression enters 
with different signs depending on whether n and m are 
both positive or both negative. This feature carries over 
to the clean case and yields 

3 

12(^) 

123 rst i=l jk 



with, 



-i'?32(^)t^93i(a=)tr iTrTln) tr (s,s]sfe) 

(2.8c) 

with C2 another positive constant. The bare values of ci 
and C2 are related, C2 = ci/16. Terms of higher order in 
q in this expansion will turn out to be irrelevant for de- 
termining the behavior at the quantum phase transition. 



C. Gaussian propagators 

We will be interested in the renormalization group 
flows of the various parameters in the field theory de- 
fined above. We will need the Gaussian propagators of 
the theory in the paramagnetic phase. These are easily 
determined from the quadratic form given by the M^, 
q^, and Mq parts of the above action. Performing a spa- 
tial Fourier transform, and using the symbol (. . .) for the 
Gaussian average, we find for the order parameter corre- 
lations 

{Arjk)M},{p)) =6k,-p5„,^mS,jlMn{k) , (2.9a) 



(r^l2(fc)i^34(P)) = -Sk-p [Sl-2-3-A- {-YSi^2A-3] 

xS,jSrsMl^2{k) , (2.9b) 

in terms of the paramagnon propagator, 

1 



Mnik) = 



f I ,,7,2 I (4Gcf/7r)|0„| 

L^uhu + \k\+GH\n„\ 



(2.9c) 



Notice that the coupling between the order parameter 
field and the fermionic degrees of freedom has produced 
the dynamical piece of M that is characteristic of clean 
itinerant ferromagnets. 

For the fermionic propagators we find. 



"rfiilfc) = '5i3^24 2?l-2(fe) 



-^i_2,3-4 27rTGKsVi-2ik) v[%{k) 



(2.10b) 



'■'^'rg34(fc) - (Sl3<524l?l-2(fe) 

-5l_2,3-48TGc2(Pi_2(fc))'Xl-2(fe) 



where V and are the propagators 

^"^^^ " iki + GHn^ 



1 



\k\ + G{H + Ks)nr. 



(2.10c) 



(2.10d) 



(2.10e) 



('2') — 1 ('2') 

Notice that "^F^ is actually the inverse of °F 



given by Eq. (|2.6b|l . while the analogous statement for 
i,2,3pi. jg j^Q^ true. This is because the coupling be- 
tween M and q gives an additional contribution to the 
fermionic spin-triplet propagator. 

Finally, due to the coupling between M and q we also 
have a mixed propagator, 

{lq,,{k)ib,,{p)) ^ [<5l-2,3-4 + (-)'^+^<5i_2,4-3] 

x,5,,%GciVri?i_2(fc)Xi-2(fc) . (2.11) 



D. Higher order terms, and diagram rules for a 
loop expansion 

The action defined by Eqs. H2.6|l - H2.8|l suffices to ex- 
tract the information we are interested in, but it is in- 
complete from a calculational point of view. Namely, in 
order to set up a loop expansion and renormalize the 
vertices in our action to one-loop order, one needs the 
term of order q^ . Although it is possible to determine 
the desired renormalizations without knowing this term 
explicitly, see below, for completeness and later reference 
we here give such a term that satisfies basic symmetry 
requirements. 

On general grounds, and by analogy with the disor- 
dered case,'^ this term must have the structure 



1,2 '■I. '■2 ii,i2 



3,4 '-3,^4 iji^ 



Xtr {Tr^TlTr.,Tl^ tr {s^^s\^S^^s\}| 

y.Y^yl{xx - X4, X2 - Xi, X3 - X4) 



x;\<zi2(^i) ^2932(^2) ;i<z34 



{X3)l\q^^{x4) 



(2.10a) 



(2.12a) 



6 



'(2,0) 



X = X 



+ 



+ 



FIG. 1: The noninteracting part of the two-point g-vertex 
to one-loop order. Solid lines denote the interacting part of 
the g-propagator, and the wave line denotes the A-propagator. 
See the text for further explanation. 



The vertex function F'^' can be expressed in terms of the 
two-point vertex r(2), Eqs. (|2.6bli>.fe|) .^ In Fourier space, 
and neglecting K^, it reads 



rlt^(fci,fc2,fc3) 



-.(2,0) 



(fci+fc2)+rif ^(fc2+fc3) . 

(2.12b) 

In addition there are terms of order and that are 
proportional to K^, as well as terms of higher order in 
but they will not be important for our purposes. 

As the last step in defining our effective field theory, 
we need to remember that setting up a g-field theory 
requires a Lagrange multiplier field A that constrains bi- 
linear products of the underlying fermion fields to the 
classical matrix field Q. In clean systems, the A-field is 
soft with a propagator that is given by minus the nonin- 
teracting part of the q-propagator,'^ 

G 

(2.13) 

This field couples to q in a way that results, upon inte- 
grating out A, in the following diagram rules'^ 

Rule 1. For calculating propagators in a loop expan- 
sion, all internal g-propagators must be taken as the in- 
teracting part of the Gaussian propagator, i.e., as the sec- 
ond term on the right-hand side of Eq. (|2.10bl) or (|2.1Uc|l . 

Rule 2. For calculating vertex functions. Rule 1 also 
applies. In addition, one needs to consider all reducible 
diagrams (which normally do not contribute to the ver- 
tices), with all reducible propagators replaced by the A- 
propagator, Eq. H2.13|l . 

As an illustration, wc show in Fig. ^ the diagrams for 
the renormalization of F^^'''^ Eq. (|2.6c|l . to one-loop or- 
der. 

This completes the definition of our effective action, 
and we will now proceed to discuss the ferromagnetic 
transition it describes. 



III. GENERALIZED MEAN-FIELD THEORY, 
AND THE FIRST-ORDER PHASE TRANSITION 

In this section we derive a generalized mean-field the- 
ory for the ferromagnetic transition in low-temperature 
itinerant electron systems. It structurally maps onto 



the generalized mean-field theory for the superconduct- 
ing transition at finite temperature.'^ The transition pre- 
dicted by these theories is of first order. We then discuss 
the conditions under which this result is stable. We will 
see that, contrary to the usual concepts concerning first- 
order phase transitions, the mean-field description can be 
invalidated by fluctuation effects that drive the transition 
second order. Physically, the first-order transition turns 
out to be unstable when it is too close to a second-order 
transition with sufhciently strong fluctuations, otherwise 
it is stable. 



A. Generalized mean-field theory 

An effective action, >loff[-M'], involving only the magne- 
tization order parameter, can be obtained by integrating 
out the fermion fields, 



„-4[JVf,g] 



(3.1) 



Here A is the action given by Eq. (|2.1bll . In general the 
evaluation of this expression is very difficult. However, 
it can be evaluated exactly within a generalized mean- 
field approximation, which is defined as follows. First, 
we ignore temporal and spatial variations of M , 



(3.2) 



Second, we restrict ourselves to Gaussian or quadratic 
order in q. That is, we replace the full action A by Eqs. 
H2.2I2.6I I2.8|l , and in these expressions we replace M by 
Eq. K^ . 

With the bare Gaussian action as given in Sec. IIIBI 
and taking into account the diagram rules mentioned at 
the end of Sec. Ill Dl there is no coupling between the soft 
modes and the order parameter. However, one needs to 
acknowledge that under renormalization, the action Aq 
will acquire a spin-triplet interaction that is generated as 
long as Ks ^ 0. Let the such generated interaction con- 
stant be Kt . Then the fermionic 2-point vertex function, 
Eq. I|2.6b|l . gets generalized to 

'rS34(fc) = <5l3^24(|fc|+Gi/|r!i_2|) 

+ '5i_3,2-4 6iQ 2ttTGKs 

(2.6b') 

This renormalization-generated spin-triplet interaction 
leads to a coupling between the soft fermionic modes and 
the order parameter in the free energy. In the resulting 
generalized mean-field approximation one obtains for the 
free energy density, /(m) = —TAcs/V, 

f = f[m^O]+t + u m'^ 

^ ^ T^lniV(A;,17„;m) , (3.3a) 



V 



k<A n 
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where A is an ultraviolet momentum cutoff, and 

+(|fc| + GHnn)^{\k\ + G{H + kt)n^f 

(3.3b) 

Minimizing / witfi respect to the magnetization gives the 
equation of state, 

h = 2tm + Aum^ 

(3.3c) 

with h an external magnetic field. 



B. Discussion of the generalized mean- field 
equation of state 

We start with some general comments regarding the 
result, Eqs. (|3.3(l . The last term in both Eqs. (|3.3a(l 
and 13.3c|) arises from fermionic fluctuations, namely, 
the Iq with r — 0,3 and i — 1,2, that are massless in 
the paramagnetic phase, but that become massive in the 
ordered phasei As discussed elsewhere,' these fluctu- 
ations lead to long-range correlations in paramagnetic 
metals, and to nonanalyticities in either the temperature 
or the wavenumber depence of correlation functions, for 
example, the magnetic susceptibility. It is also interesting 
to note that Eqs. 1)3. 3|) are identical to the equations de- 
scribing the first-order phase transition in conventional 
superconductors at finite temperature.' As mentioned 
in the Introduction, the physics of the respective phase 
transitions is very similar as well. 

With some work, the integrals or sums in Eqs. (|3.3|) 
can be explicitly performed. However, the most impor- 
tant features can be obtained by inspection and simple 
asymptotic analysis. At zero temperature, and for small 
m, the leading nonanalytic m-dependence is a negative 
term on the right-hand side of Eq. (|3.3c|l that is of or- 
der m'^ in generic dimensions, and of order In 1/m in 
d = 3. At low but finite temperatures this nonanalyticity 
is effectively replaced by a negative term of order 

m(m2 + const, x T'^)('i'i)/2 

in generic dimensions, or 

TOMnl/(m2 + const, x t2)1/2 

in d = 3.^ Here const, is a positive constant proportional 
to H^{1 + H / KtY / c^. Analogous terms, with an extra 
factor of TO, appear in Eq. (|3.3a|l . As T — > Eq. (|3.3a(l 
thus has the standard form of a free energy, or effective 
potential, that leads to a discontinuous phase transition 
at some t = ti > Q, see Fig. |21 Schematically, this free 



f 




m 



FIG. 2: Schematic form of the free energy as a function of 
the order parameter. 




t 

FIG. 3: Schematic form of the pheise diagram at ^ = 0. 
The dashed line denotes a second-order transition, the solid 
line denotes a first-order transition, and TCP denotes the 
tricritical point. 

energy functional in the presence of an external magnetic 
field h can be written, in 1 < d < 3, 

/ = /(TO = 0) + im2-i;TO2(m2 + r2)(''-l)/2 

-fum — hm + . . . (3.4a) 

and in d = 3, 

/ = f{m = 0) + tm^ + vm^\n{m'^ +T'^) 

+um'^ - hm + . . . (3.4b) 

In this schematic representation, the mean-field equation 
of state in the most interesting case, d — i, takes the 
form 

h 2^TO + 4^;m^ln(TO^ -l-T^) 

+ -^(4- + 2.-^) . 

(3.4c) 

In these equations we use units such that /, to, and T 
are measured in terms of a microscopic energy, e.g., the 
Fermi energy, t, v, and u are then all dimensionless. 
w > is quadratic in Kt or cf, so in strongly correlated 
systems v is larger than in weakly correlated ones. 

In d = 3, these equations predict the phase diagram 
shown in Fig. |31 There is a tricritical point at 

T = Ttc = exp(-u/2i;) . (3.5a) 
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At T = 0, there is a first-order phase transition at t = ti, 
with the magnetization changing discontinuously from 
zero to a value mi. One finds 

nil = exp[-i(l + u/u)] , (3.5b) 
ti — vml . (3.5c) 

In d = 2, there is no finite-temperature magnetic phase 
transition. However, at zero temperature there is a QPT, 
which is predicted by the Eqs. (|3.4|l to be discontinuous. 
The discontinuity in the magnetization and the transition 
point are given by 

mi = {3v/Auf , (3.6a) 
ti = ^m^ . (3.6b) 

In d > 3 the nonanalyticitic terms produced by the soft 
modes are subleading, and the transition is described by 
ordinary mean-field theory. The generalized mean-field 
theory thus suggests an upper critical dimension = 3. 
As we will see in the next section, a more sophisticated 
analysis confirms this result. 

C. Validity of the mean- field description 

Normally, first-order phase transitions are not sensi- 
tive to fluctuation effects. We now argue, however, that 
in the present case of a first order transition driven by 
fiuctuations (viz., soft fcrmion modes), additional fluc- 
tuation effects can destabilize the mechanism underlying 
the flrst-order transition and lead to a fluctuation-driven 
second-order transition. This will happen if, in a sense 
described below, the first-order transition is too close to 
an unrealized second-order one. To illustrate this point, 
consider the two schematic free energy functionals shown 
in Fig. 21 In Fig.0Ja) the first-order transition occurs far 
from t = 0, and fluctuation effects are negligible. How- 
ever, for the case shown in Fig.QJb) the fluctuations near 
the (unrealized) second-order transition at t = can af- 
fect the first-order transition that preempts the second- 
order one, and need to be taken into account. 

Before presenting technical details in the next section, 
let us elaborate on this general point. First, we note 
that as long as one is far from any continuous transi- 
tion (which brings in new fiuctuation effects) the func- 
tional forms of the free energy funtions given by Eqs. 
()3.4(l are exact for small magnetizations. This follows 
from the properties of a Fermi-liquid fixed point and the 
corrections to scaling near it.'^ The mean-field descrip- 
tion above suggest a second-order, or continuous, phase 
transition at t — 0, which is preempted by the first-order 
transition. The latter results from fluctuations that are 
germane to a Fermi liquid and have nothing to do with 
any critical point. These observations imply that the 
first-order transition discussed above will take necessar- 
ily take place if (1) toi is small enough, and (2) ti is 



f 




m 



FIG. 4: Schematic forms of the free energy as a function of 
the order parameter. In (a) the first-order transition is not 
affected by fluctuations, in (b) it may be. 



not too small, so that additional fluctuation effects due 
to the underlying critical point at t = can be ignored. 
Examining the Eqs. (|3.4|) and H3.5|l we see that this can 
occur when both v (which describes correlation effects) 
and u/v are large. More generally, it is reasonable to 
expect a first-order phase transition, with no restrictions 
on m, whenever correlation effects are large. 

Second, the next natural question is, what happens if 
this is not the case? In particular, we note the following. 
The Eqs. H3.3|l imply that the coefficient v is inversely 
proportional to H. As explained in Section Hvl below. H 
is proportional to the specific heat coefficient. Since H is 
expected to be sensitive to critical fiuctuations, and per- 
turbation theory suggests a logarithmic divergence at a 
continuous transition in d = 3, this suggests that v might 
be suppressed close to, and at, a continuous transition. 
To examine this possibility one needs to go beyond simple 
perturbation theory. In the next Section renormalization 
group methods are used to untangle the possibilities. We 
find that a continuous transition does indeed occur if ti 
is sufficiently small, and if the bare u is sufficiently large 
to stabilize the effects of critical fermionic fluctuations 
that are otherwise suppressed. 



IV. RENORMALIZATION GROUP FLOWS, 
AND THE SECOND-ORDER PHASE 
TRANSITION 

A. Renormalization of the effective action 

The parameters t, a, u, G, H, ci, C2 as well as the fields 
M and q in the theory defined above are renormalized 
under renormalization group (RG) transformations. We 
will employ a differential momentum-shell RG and inte- 
grate over all frequencies. If b is the RG length rescaling 
factor, then we rescale wavenumbers and the two fields 
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k ^ bk , 



(4.1a) 
(4.1b) 
(4.1c) 



Here rjq and rjM are exponents that characterize the 
spatial correlations of the fermion and order-parameter 
fields, respectively. The rescaling of imaginary time, fre- 
quency, or temperature is less straightforward. We need 
to acknowledge the fact that there are two different time 
scales in the problem, namely, one that is associated with 
the critical order-parameter fluctuations, and one that is 
associated with the soft fermionic fluctuations. Accord- 
ingly, we must allow for two different dynamical expo- 
nents, zm and Zq, and imaginary time and temperature 
may get rescaled either via 



T b'"-'T 



or via 



T^b-'-'^T , T-^b^'-'T , 



(4.1d) 



(4.1e) 



How these various exponents should be chosen is dis- 
cussed below iS 



1. Zero-loop flows 

In the tree, or zero-loop, approximation the RG flow 
equations for the parameters in our field theory are easily 
determined by power counting from the action given in 
Sec. iHBl With ^ = In 6 we find 



da 
du 



(4.2a) 
(4.2b) 
(4.2c) 
(4.2d) 
(4.2e) 
(4.2f) 



-77 = -{d + ZM + 2ryM - 4) M , 
ai 

dci 1 

^ - \{&-d-z-2i^q-i^M)c2 . (4.2g) 

Now we examine these flow equations in order to see 
whether the allow for a critical (i.e., unstable in only one 
direction) fixed point, at least above some upper critical 
dimension. This will amount to an analysis of the sta- 
bility, or otherwise, of Hertz's critical fixed point. Note 
that in giving Eqs. Ij4.2f|l and | |4.2g| |, the particular choice 
of z was not yet specified because it is not obvious if a Zq 
or a zm should be used for these terms that describe a 




FIG. 5: Example of a fermionic loop renormalizing the para- 
magnon propagator, with the dashed line representing the 
magnetization field. Since the loop integral is over fermionic 
propagators only, the C2 associated with the vertices carry a 
time scale given by Zq. 



coupling between q and M fields. For the analogous dis- 
ordered electron problem, this point has been discussed 
in detail in Ref. ? . 

If we assume the Fermi-liquid degrees of freedom to be 
at a stable Fermi-liquid fixed point, we must choose G 
and H to be marginal, which implies 



1 



1 



(4.3) 



Hertz'^ further chose (in our language) a and ci to be 
marginal, which implies 

7?M = , (4.4a) 

and 

ZM = 3 . (4.4b) 

The latter choice is motivated by the paramagnon prop- 
agator, Eq. H2.9c|) . which yields appropriate mean-field 
critical behavior only for a marginal ci, given that G 
and H are marginal. This also implies that 2; = zm in 

Eq. gisgi. 

With these choices, t is the relevant variable charac- 
terizing the critical fixed point, and Eq. H4.2a|l yields a 
correlation length exponent v = 1/2. The variable u is 
irrelevant for d > 1, suggesting an upper critical dimen- 
sion equal to unity. Indeed, Hertz's conclusion was that 
the mean-field fixed point characterized by the above ex- 
ponent was stable for c? > 1. 

However, we still need to examine the behavior of C2. 
It is irrelevant for c? > 1 if we use z = zm in Eq. | |4.2g| ). 
However, as already emphasized in Ref. ? , one also 
has to consider the case z = Zq 'm this equation. This 
becomes obvious if one uses the M-q^ vertex, whose cou- 
pling constant is C2, to construct loops. Clearly, pure 
fermion loops appear, the simplest example of which is 
shown in Fig. [S] and in this case z — Zq is the appropri- 
ate choice. We illustrate this point below by means of an 
explicit calculation. 

Using z — Zq in Eq. ( |4.2g| ), we see that C2 becomes 
relevant with respect to Hertz's fixed point for d < 3, 
giving an actual upper critical dimension = 3. This 
is in agreement with the result of the generalized mean- 
field theory, see Sec. IIII Bl Physically, this surprising 
results means that soft or gapless fermion excitations play 
an important role in determining the phase transition 
behavior at, and below, three spatial dimensions even 
though naive power counting suggests d+ = 1. This is 
further examined in the next subsection, as well as in 
Section HVBl 
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6a = 



5u = 



6H = 



Q 



FIG. 6: Diagrams that give the leading renormahzation of a, 
u, and H, to one-loop order in d = 3. The quartic vertex in 
the diagram for 5H was defined in Fig. Q 



2. One-loop flows 

In this subsection we examine the effects of C2 on the 
phase transition. We will be mainly concerned with the 
behavior in three-dimensions, the behavior in d < 3 will 
be discussed using other techniques in Sec . IIV 51 below. 

In d = 3 the relevant diagrams can, in principle, give 
logarithmic corrections or renormalizations to the various 
coupling constants. Taking into account that there are 
two time scales, it is easy to show by power counting that 
there will be no logarithmic corrections to ci, C2, G, and 
t. This implies that for these coupling constants, the 
flow equations given in Eqs. 1)4.2(1 remain valid to one- 
loop order. Motivated by the disordered case, we will 
be looking for a fixed point where G is marginal, which 
implies 



Vq = 1 



(4.5a) 



We further require ci (with z — zm) and C2 (with z — Zq) 
to be marginal, which implies 



and 



Zm + Vm 



Zq + riM = 1 



(4.5b) 



(4.5c) 



Of the various scale dimensions introduced above, this 
leaves only one, e.g. t^Mi as independent. For the irrele- 
vant version of C2 (with z = zm), Eqs. H4.5a|l and (|4.5b|l 
imply 



dc2 



-C2 



(4.6) 



For the remaining quantities, power counting shows 
that they do allow for logarithmic renormalizations in 
d = 3. The diagrams that give rise to these are shown 
in Fig. |B1 A crucial feature is that u is renormalized 
by a negative logarithmic term. In a purely perturbative 
treatment, this implies that u changes sign, which in turn 
implies a fluctuation driven first-order phase transition. 



and the existence of a tricritical point at finite tempera- 
tures, consistent with the generalized mean-field theory. 
However, the renormahzation group flow equations re- 
sum perturbation theory in a specific way, and in this 
subsection we show that this 'tricritical' behavior does 
not necessarily persist to all orders. 

The explicit flow equations are obtained by evaluating 
the diagrams shown in Fig. Determining the general 
structure of the flow equations does not require a detailed 
calculation, but can be achieved by power counting. At 
zero temperature, we find 



da 
M 
du 

di 
dH_ 

m 



-TjMa 



H 



= -{2 + riM)u-Au 



H 



VmH 



A 



H 



a + t 



(4.7a) 
(4.7b) 
(4.7c) 



where the At are positive constants. In giving Eqs. 1(4.7(1 
we have absorbed the marginal coupling constant G and 
the marginal version of C2 into these constantSjiS 

The prefactors Ai can be determined by a detailed cal- 
culation of the diagrams. In the case of Aa and Ah, one 
can also obtain the result by the following alternative 
method, a is the coefhcient of the gradient squared term 
in the spin susceptibility of a non-magnetic reference sys- 
tem. The logarithmic renormahzation of the latter in 
d = 3 has been calculated in Ref. ? , and we can thus 
find the renormahzation of a from that paper. For H, 
we notice that it is related to the specific heat coefficient 
Jv = Gv/T by 



7y = 8ttH/3 



(4.8) 



This relation between the frequency coupling constant 
and the specific heat was first established for disordered 
electron systems by Castellani and Di Castro.'^ A proof 
by means of Ward identities'^ applies to clean systems 
as well. One can therefore obtain the renormahzation of 
H from a calculation of the specific heat, which in turn 
follows from the Gaussian free energy density Jq. From 
Sec. Ill Cl we find, at criticality, 



3r 



/G = i7EEl-(«P'+4Gc?AIPl) 



(4.9) 



and the specific heat coefficient is obtained by differenti- 
ating twice with respect to temperature. 

We have chosen the second method to calculate Aa and 
Ah, and obtair 



,11 



Aa = GcI/9tt^ 
Ah = SGcI/tt^ 



(4.10a) 
(4.10b) 



For later reference we note that An/Aa — 27 > 1, but 
for now we consider the general case. The value of A^ 
will not be needed, other than that it is positive. 
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Next we solve the Eqs. (|4.7ll . first for the case t = 0. 
It is convenient to first construct flow equations for the 
quantities 

f^cl/H , g^aH . (4.11a) 
The flow equations for these objects in d = 3 are, 
dg 



dl 
dl 



= Ah - Aa , 
= -{2 + mi)f -Anf/g 



(4.11b) 
(4.11c) 



Solving these equations and using the result in Eq. (|4.7b|l 
gives 



with 

K^2 + riM, A^ 
and 



Aufo 
A{k-1) 



Ah - Ag 
90 



1 



1 



(i + Ae) 



K 



K-l 



(4.12a) 



A 



H 



90 = g{i = 0) = l/96^«F 

fo = f{e = 0) = 



Ah - Aa 

(4.12b) 



(4.12c) 
(4.12d) 



Since K > 1 and A > for Ah > Aa, we see that in 
contrast to the perturbative result, u{£) does not neces- 
sarily become negative as ^ — > (X). Rather, the term in 
braces in Eq. H4.12a|l asymptotically approaches a value 
Mo — Aufogo/ Aa- Depending on microscopic parameter 
values, u thus may or may not become negative for large 
scales. We conclude that a nontrivial continuous phase 
transition may exist for d < 3. 

We note, though, that if the opposite inequality were 
to hold. Ah < Aa, than A would be negative and u 
would become complex at large scales. These two fea- 
tures would indicate a first-order phase transition. This 
suggests that the actual first-order phase transition oc- 
curs at i > where Ah is effectively smaller, cf. Eq. 



B. Critical behavior 

1. Critical behavior in d — 3 

We are now in a position to determine the critical be- 
havior at the second order transition that we have shown 
in the previous subsection to exist in a certain regime of 
parameter values. In d = 3 we do so by using the explicit 
solution of the flow equations given in Sec. IIV A 21 above. 

Let us consider the paramagnon propagator in the crit- 
ical regime. Since t, c\, and G are not singularly renor- 
malized at one-loop order, while the coefficient a acquires 



a nontrivial renormalization, it has the form, see Eq. 
(EScJ, 



M.(fc) = l/(t + a(fc)fc2-|-|f^„|/|fc|) 



(4.13a) 



fe and ri„ have been made dimensionless by means of suit- 
able microscopic scales. The fc dependence of a follows 
from Eqs. (|4.t)|l and (|4.11|) once r\M has been chosenii^ In 
writing Eq. (|4.13a|l we have tacitly assumed that there is 
no 'wavefunction renormalization' that would give the 
numerator a scale dependence. Since r\M determines 
the scale dimension of the magnetization field, see Eq. 
(j4.1bp . and hence the wavefunction renormalization, we 
need to choose r\M = in order to be consistent with 
this assumption. From Eqs. (|4.6() and H4.11|l in the limit 
£ ~ In l/|k| — > cx) we then obtain 



a(fc ^ 0) cx (lnl/|fc|) 



-1/26 



(4.13b) 



Such logarithmic corrections to power-law scaling can be 
conveniently expressed in terms of scale dependent criti- 
cal exponents. For instance, with 6 l/|fc| a RG length 
scale factor^^ we can write a(fc) oc |fcp^'', with a scale 
dependent critical exponent 77 given by 



77 = — In In 6/ In 6 
' 26 ' 



(4.14a) 



We stress that 77 is the physical critical exponent that de- 
scribes the wavenumber dependence of the paramagnon 
propagator at criticality, M. oc |fc|~^+'', as opposed to 
r\M, which has no direct physical meaning. In Appendix 
^we demonstrate that a different choice of r\M leads to 
the same physical result. 

The correlation length exponent v, the susceptibility 
exponent 7, and the dynamical exponent z can be di- 
rectly read off Eqs. H4.13|l . viz. 

v^\l(2-T]) , z^3-r] , 7=1 . (4.14b) 

These exponents are defined as usual, i.e., ^ cx t^'^, 
il T ^ A4 oc , with ^ the correlation length. 
The physical dynamical exponent z is different from the 
exponent zm in Eq. 



Ij4.5b|l , for the same reason for which 
V 7^ Vm- Also notice that l/i/ is not given by the scale 
dependence of t that results from the i-flow equation, Eq. 
(|4.2all . since the scale dependent coefhcient a is a dan- 
gerous irrelevant variable with respect to the correlation 
length. 

The order parameter exponents /3 and S (defined by 
m (X and 771, oc h-^^^, respectively, with h a magnetic 
field) can be obtained from scaling arguments for the free 
energy, see Appendix We find 



P=l/2 



(5 = 3 



(4.14c) 



Finally, we define a specific heat exponent a by Cy oc 
T"" at criticality^ It can be determined by either of 
three methods, viz., (1) Eq. (|4.8|) together with the solu- 
tion of the flow equation for H, (2) renormalized pertur- 
bation theory for the free energy, i.e., Eq. (j4.13b|) in Eq. 
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(|4.9|) . or (3) a scaling argument for the free energy, see 
Appendix^ Either way we obtain the exact relation 

a = -1 + (lnln6/lnfe- ?7)/z . (4.14d) 

The result for ry is valid to leading logarithmic accu- 
racy; the values of 7, (3, and 5, as well as the relations 
between 77 and v, z, and a, respectively, are exact. 

2. Critical behavior in d < 3 

In dimensions less than three, the critical behavior can 
be controlled by means of an expansion in e = 3 — d. We 
are again looking for a fixed point where G and ci are 
marginal, so Eqs. H4.5a|l and Ij4.5b|l still hold. Equation 
(I4.5cll gets generalized to 

Zq + T]M^l + e , (4.15) 

which guarantees that C2 with z — Zq is still marginal. 
We then look for a fixed point where a and H are both 
marginal. rjM then coincides with the physical exponent 
rj, as it does in the alternative treatment of the case d = 3 
given in Appendix^ We find 

77 = -e/(AHM„-l) = -6/26 + 0(e2) . (4.16a) 

The other exponent follow from this, i^, z, 7, /3, and 6 
are still given by Eqs. l)4.14bl4.14c|l . and for the specific 
heat exponent a we have 

a = -d/(3-77) . (4.16b) 

V. SUMMARY AND DISCUSSION 

We summarize the achievements of this paper as fol- 
lows. 

First, we have given an effective field theory that de- 
scribes the quantum ferromagnetic transition in clean 
electronic systems. It involves coupled fields describing 
the magnetization degrees of freedom, as well as gapless 
fermionic excitations. If the effects of the latter are ne- 
glected beyond tree level, as was the case in earlier theo- 
ries describing this quantum phase transition,' then the 
resulting description of the phase transition is incorrect 
for all d < 3. That is, the coupling to the fermionic de- 
grees of freedom leads to an upper critical dimension for 
this phase transition of d+ = 3. 

Second, we have shown that the fermionic fluctua- 
tions lead to two very different types of fluctuation-driven 
quantum phase transitions, depending on microscopic de- 
tails. Generically, the quantum ferromagnetic transition 
in d < 3 is a fluctuation-driven first-order transition. 
This is in contrast to the conventional result,'^ as well 
as to the Landau theory description of this phase tran- 
sition. In d = 3 we have also discussed the situation at 
low, but finite, temperatures. In general we argue that 



this system will have a tricritical point separating lines 
of second and first-order phase transitions. A schematic 
phase diagram is shown in Fig. 13 These results are in 
agreement with the experimental observations in MnSi'^ 
and UGe2.^ 

Third, we have shown that if the microscopic details 
are such that the fluctuation-driven first-order quantum 
phase transition is too close to a second-order, or contin- 
uous, transition, then critical fluctuations will suppress 
the fermionic fluctuation effects that lead to a first-order 
transition, and a fluctuation-driven second-order transi- 
tion resultsi^ For this case, the critical behavior in d = 3 
has been computed and has been found to be mean field- 
like, with logarithmic corrections. For d < 3, the critical 
behavior is nontrivial, but can be controlled by means of 
a 3 — e expansion. Both the possibility of a second-order 
transition and the fact that the critical behavior in this 
case is essentially mean field-like is in agreement with the 
experimental observations on TxLw-i and NixPdi_x.'^ 
Our theory thus explains a rather confusing experimen- 
tal situation, where the transition in bulk systems is ob- 
served to be continuous in some systems, and discontinu- 
ous in others. We further note that the fluctuation effect 
that leads to a first-order phase transition grows with 
the strength of electronic correlation, or interaction, ef- 
fects. This suggests that in strongly correlated systems 
a first-order transition is generally expected. 

Fourth, we have noted a mathematical and physical 
relation between the fluctuation-driven first-order phase 
transition discussed here, and the ones known to occur in 
finite-temperature superconductors and in liquid crystal 
systems.'^ ^' In all these systems, soft modes couple to 
the order parameter fluctuations in such a way that their 
contribution to the free energy is reduced in the ordered 
phase. It is this mechanism that causes the discontinuous 
transition to occur. The fluctuation-driven second-order 
transition discussed here is similar to the one that occurs 
in classical Potts models.'^ 

Elsewhere we have discussed the effects of nonmag- 
netic disorder on this phase transition, and on the phase 
diagram shown in Fig.|21'^ In general, sufficiently strong 
disorder drives the tricritical point shown in Fig.|21to zero 
temperature, making the zero-temperature transition in 
the presence of sufficiently strong disorder continuous. 
This quantum phase transition is in a different universal- 
ity class than the fluctuation-driven second-order tran- 
sition in clean systems discussed above, and its critical 
behavior has been determined exactly.' 
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APPENDIX A: RG WITH SCALE DEPENDENT 
EXPONENTS 



Our choice of r/M = in Sec. IIV B II is somewhat un- 
conventional since it makes the coefficient a irrelevant, 
rather than marginal. In this Appendix we demonstrate 
that making a marginal leads to the same physical re- 
sults, but comes with complications of its own. 

Since the physical exponent 77 is scale dependent, see 
Eq. (|4.14a() . choosing a marginal requires a scale depen- 
dent r/M- This changes the flow equations. If we still 
require that G, ci (with z = zm) and C2 (with z — Zg) 
are marginal, we obtain, instead of Eqs. (|4.2II4.7I1 . 



da 

di 
dH 

m 

dc2 

m 

du 

dl 



f]M{l)a- Aa/H 
Ah 



a + 1 



-C2 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
(Ale) 



Here C2 in Eqs. IjAldp and (|Ale|l refers to the irrelevant 
(z = Zm) incarnation of C2, and 



drj 



(A2) 



We now look for a fixed point with respect to which C2 
and u are irrelevant, while a is marginal. The latter 
condition leads to a differential equation for rjM, 



dr]M 
dl 



1 - K 



= 



with K from Eq. H4.12b|l . In deriving Eq. HA3|) we have 
used the fact that g — aH obeys a fiow equation, aXt — 0, 



dg 



An-Aa 



(A4) 



Since a is marginal, 77M now represents the physical expo- 
nent 77. Eq. (|A3I) is easy to solve, and for asymptotically 
large values of i one recovers Eq. (|4.14a|l . Similarly, the 
correlation length exponent v is now given by the scale 
dimension of the relevant operator i, and from Eq. (|Ala|) 
we recover the first equality in Eq. (|4.14b(l . It is easy to 
check that all other physical results also agree with Sec. 
IIV B II 



APPENDIX B: SCALING ARGUMENTS FOR 
THE FREE ENERGY 

In this Appendix we consider the scaling behavior of 
the free energy density. Let us add a magnetic field term 



Ah our action. 



Ah 



h - dx M{x) 



(Bl) 



h gets rescaled via h b^'^^h with a scale dimension 



[h] = - {d + 2 + ZM -Vm) 



(B2) 



and the free energy density obeys the scaling law 

f{t,T,h) = h-^'^+'''^ f{th^''',Th''''' ,hb^^^) . (B3) 

The magnetic susceptibility is given by x = f /dh?, 
and it is readily checked that Eq. l|B3p reproduces the 
susceptibility exponent 7=1 that we obtained in Sec. 
IrVBl from the paramagnon propagator. The critical be- 
havior of the magnetization, m = df /dh can be obtained 
from Eq. ljB3|l as well if we take into account that the ir- 
relevant variable u is a dangerous irrelevant operator with 
respect to m (but not with respect to %)■ We thus need 
to include u in the set of scaling variables, and find 



i(t, h, u) = 6-('i+^A/- m (t b^/'',h , u 



with 



-{d-l + riM) , 



(B4a) 



(B4b) 



the scale dimension of u, see Eqs. (I4.2cl I4.5c|l . Taking 



(A3) 

into account m{h = 0) oc u ' and m{t ~ 0) 



we obtain mean- field values, Eq. (|4.14c|l . for the order 
parameter exponents P and 6. 

We also comment on the relation between the free en- 
ergy scaling and the specific heat exponent a. From Eq. 
(|B3|) we obtain a scaling law for the specific heat coeffi- 
cient ■^v — d^f/dT^ at criticahty, 



IV {T) 



(B5) 



In d = 3 — e this agrees with the result for the expo- 
nent a, Eq. (|4.16b|) . as obtained from either the H-How 
equation, or renormalized perturbation theory. In d = 3, 
however, the simple scaling argument does not agree with 
the other two methods. The reason lies in the fact that in 
d = 3, d — Zm apart from logarithmic corrections to scal- 
ing. This is one of the 'resonances' between exponents 
that have been discussed by Wegner in his classification 
of sources of logarithmic corrections to scaling.'^ This 
resonance leads to an additional logarithm that is missed 
by the simple scaling argument. Once this is taken into 
account, all three methods agree in d = 3 as well. 
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By "clean" we mean a system that is free of impurities or 
structural disorder. 

The sign of this term is consistent with the fact that the 
mode-mode coupling contribution is due to fluctuation ef- 
fects that tend to suppress ferromagnetism. In a system 
with quenched disorder, the effect has the opposite sign. 
See Ref. ? for a discussion of this point. 
Of course the magnetization is also if electronic origin. 
However, since the order parameter fluctuations play a spe- 
cial role in the theory, we separate them out and refer to 
all other soft modes as "fermionic". 

Both of these approaches have been widely used in the 
history of phase transitions. For instance, the Ginzburg- 
Landau theory of superconductivity was originally deduced 
from general symmetry arguments and only later derived 
from the microscopic BCS theory. 

We deliberately use the same notation as in our treatment 
of disordered ferromagnets, Ref. ? , to underscore the 
similar structures of the two theories. 

This structure can be deduced as follows. Integrating out 
the massive modes in tree approximation proceeds formally 
as in the disordered case, where it produces a nonlinear 
sigma model, see Refs. ? ? . The only difference is that 
in the clean case, soft single-particle excitations have been 
integrated out to produce the effective Q-field theory. This 
leads to a singular vertex in the gradient-squared term of 
the sigma model that is, in the long-wavelength limit, pro- 
portional to an inverse wavenumber and changes the dis- 
persion relation of the soft modes from a diffusive one to 
a linear one, see Eq. 12.6cl . The requirement that the two- 
point vertex F'^' remains soft under renormalization puts 
constraints on the higher vertices in the g-expansion. These 
constraints are fulfilled if the singular vertex is the same 
for all terms in the g-expansion. This suggests that the 
vertex, if expressed in terms of the vertex, is the same as 
in the nonlinear sigma model. This procedure leads to Eqs. 
(ITT^ . In Sec.|IV|we wiU see that the results obtained from 
this F'*-* agree with those obtained by other, more indirect, 
means. We also mention that the above arguments suggest 
that one can construct an effective field theory for the soft 
modes in a clean fermion system that is analogous, and 
closely related, to the well-known nonlinear sigma model 
that describes disordered fermions. This general theory will 
be pursued separately. 

This can be seen explicitly by using Eq. 13.21 in Eq. (12. Bet 

and re-calculating the g-propagator. Notice that the mean 



magnetization m acts like an external magnetic field, which 
breaks the symmetry in spin space and gives two of the 
three soft spin-triplet modes a mass. This is the clean ana- 
log of the "spin diffusons" in the disordered case, which 
also acquire a mass in either an external magnetic field, 
or in a phase with a nonvanishing magnetization, see, e.g., 
Refs. ? ? . 

* Notice that there is no term proportional to mT^luT 
in d = 3. In other words, there is no renormalization of 
i, or the magnetic susceptibility, that is proportional to 
In T. This feature of the generalized mean-field theory 
is in agreement with exact perturbative calculations, Ref. 
? , as well as with Landau Fermi-liquid theory, Ref. ? . 

^ We use Ma's method for identifying simple RG fixed 
points. Accordingly, we use physical arguments to choose 
the values of various exponents, and then check self- 
consistently whether these choices indeed lead to appro- 
priate fixed points. See Ref. ? . 
^'^ In the case of Eq. 114. Tel . this requires some explanation, 
as the C2 that appears here is nominally the irrelevant C2 
(with 2 = zm). However, it effectively acts like a marginal 
operator since the vertex function F*-^-* is proportional to 
a frequency rather than being a constant. This mechanism 
for a nominally irrelevant operator turning into a marginal 
one is the same as in the disordered case and has been 
discussed in detail in Ref. ? . 

Here we use the fact that the marginal version of C2 is 
related to the marginal operator ci, see the remark after 
Eq. lITHct . 

We remind the reader that the choice of r]M is in principle 
arbitrary, although some choices make the analysis of the 
critical behavior easier than others. 

We use the notation o ~ 6 for "scales like 6" , and a oc b for 
"a is proportional to 6" . 

This is a generalization of the usual definition of a at ther- 
mal phase transitions. 

This second-order transition is unrelated to the one dis- 
cussed elsewhere,' and the two transitions belong to two 
different universality classes. While the one discussed in 
Ref. ? could be realized somewhere, it is not consistent 
with low-order perturbation theory, and its realization re- 
quires something qualitatively to change at higher order. 
The second-order transition discussed in the present pa- 
per, on the other hand, is consistent with everything that 
is known. 



